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Abstract. The rotational influence on the eddy-diffusivity tensor Dij for anisotropic turbulence fields is considered in order 
to explain the lithium decay law during the spin-down process of solar-type stars. Rotation proves to be highly effective in the 
transfer of chemicals through the solar tachocline (beneath the convection zone) which is assumed to contain only turbulence 
with horizontal motions. The effect is so strong that the tachocline turbulence must not exceed a limit of ~10 - ' 3 " 4 ' of the 
rms velocity in the convection zone in order to let the lithium survive after Gigayears. Such long depletion times can also be 
explained by a very small rotational influence upon the eddy-diffusion tensor if it is realized with correlation times shorter 
than 15 min. It is argued that such slow and/or short-living turbulence beneath the convection zone could hardly drive the solar 
dynamo. 

In our theory the diffusion remains small for rapid rotation due to the rotational quenching of the turbulence. In young stellar 
clusters like Pleiades, there should be a (positive) correlation between rotation rate and lithium abundance, where the fastest 
stars should have maximal lithium. First inspections of the data seem to confirm this finding. 
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! 1. Introduction 

According to Duncan (1981), the lithium at the surface of 
cool main-sequence stars decays exponentially, unlike the stel- 
■ lar spin-down, which follows a power law i -1 / 2 (Skumanich 
1972). The characteristic decay time is about 1 Gyr. The pri- 
mordial lithium is destroyed by nuclear reactions at temper- 
atures in excess of 2.6 • 10 6 K (see Michaud & Charbonneau 
1991; Ahrensetal. 1992), i.e. about 40.000 km below the lower 
edge of the convection zone (Fig. [l]). There must be a drift pro- 
cess for the chemicals from the bottom of the convection zone 
through the solar 'tachocline' to the burning domain. The ef- 
fect must be small, however, in order to allow the existence of 
lithium in the solar atmosphere even after 4.6 Gyr. The lithium 
decay time is about 10 7 times the convection zone diffusion 
time of ~ 100 yrs. We are thus looking for a rather small ef- 
fect which, however, cannot simply be microscopic diffusion 
(Schatzman & Maeder 1981; Spruit 1984; Zahn 1989, 1992). 
In the following, the consequences are presented of a quasi- 
linear mean-field approximation of an anisotropic turbulence 
field which might be located in the solar 'tachocline', which 
is stably stratified unlike the unstably-stratified solar convec- 
tion zone. However, to form reduction factors of 10 7 "" 8 for any 
eddy diffusivity is a nontrivial problem. In other words, the de- 
Send offprint requests to: G. Riidiger 




AGE [Gyr] 



Fig. 1. The outer structure of the Sun during its main-sequence 
life. Ri limits the unstably stratified outer convection zone 
while Rq displays the lithium burning zone with a tempera- 
ture of 2.6 Mio K. Computed by S. Kammerer (2000) with the 
solar model of Stix & Skaley (1990) 



sired eddy diffusion coefficient only exceeds the molecular one 
by two orders of magnitude (Baglin etal. 1985; Lebreton & 
Maeder 1987) and this is not easy to explain (Vincent etal. 
1996). 
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The transport of a passive scalar is governed by the diffu- 
sion equation, 



d P C 
dt 



div {pCu - pDVC) = 0, 



(1) 



with D as the microscopic diffusion coefficient. In the sense of 
the anelastic approximation we shall always apply the source- 
free condition of the mass flux, i.e. 



divpit = 0. 



(2) 



If the field u describes a turbulent flow pattern, the fields are 
split into a mean and a fluctuating part, 



c = c + c, 



u = u + u 



(3) 



Averaging ([I]) we get the well-known diffusion equation in the 
presence of turbulence 



d P C 
~dt 



div {p(C'u') - pDVC) = 0. 



(4) 



The influence of a meridional circulation considered by 
Charbonneau & Michaud (1991), Zahn (1992) and Chaboyer 
& Zahn (1992) is neglected here. 

The fluctuations of the chemical concentration in a quasi- 
linear approximation follow the relation 



dt 



div (pCu 1 - pDVC) = 0. 



(5) 



From this we have to compute the turbulent concentration flux 
vector (C'u 1 ) which must be inserted into (0) in order to find 
the mean-field diffusion equation. In particular, we are inter- 
ested in the effect of the basic stellar rotation and its conse- 
quences during the stellar spin-down process (see Charbonnel 
etal. 1992). 

2. The diffusion tensor 

In a corotating frame of reference the equation for the fluctuat- 
ing part of the velocity field u' is 



du' 1 

— + 2fl x u' + ~W - f t Am' = /' 

dt p 



(6) 



with /' as the random turbulence force and v t as some back- 
ground eddy viscosity. The desired correlation between con- 
centration fluctuations and velocity fluctuations, (C'u 1 ), can 
be found by Fourier transforming (ra) using 



u'(x,t) = J u{k,Lo)e i ^^dkduj. 



(7) 



As usual, the original anisotropic turbulence is assumed to be 
stationary, homogeneous and anisotropic, i.e. 

(uf (k, cj)uf{k', «')> = <K W + ^') S ( k + k ') 
[ E{k,u) , s 3 
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E 2 (k,iu)[5 lJ -k°k° - 



-(g-k°) Sij - g % gj + {g ■ k°){g t k° + g k°) 



(cf. Riidiger 1989). E(k,u>) is the isotropic part of the turbu- 
lence spectrum, fc° = k/k and g is the vertical unit vector. The 
radial turbulence intensity may be denoted by w 2 — (u^ 2 ) 
while for the azimuthal turbulence intensity v 2 = {u^ 2 ). An 
anisotropy parameter s is defined by v 2 — sw 2 , so that a large 
s denotes horizontal-type turbulences. It is 



E(k, to)dkduj. 



(9) 



The vertical vector g represents the basic anisotropy, which is 
described by the spectrum E^ik, oj) of additive horizontal mo- 
tions 



v 2 = - J J E 2 {k,uj)dkduj. 
o 



(10) 



The result for the turbulent concentration-flux vector may be 
written as an anisotropic diffusion in terms of the mean con- 
centration gradient, i.e. 



(cf. Dolginov & Silantev 1992). We shall compute in the fol- 
lowing the diffusion tensor without rotation and with rotation 
for the isotropic and anisotropic parts of the turbulence fields. 

2. 1. Isotropic turbulence 

The result is very simple without rotation. It follows 



D; 



1 



Dk 2 



D 2 k 4 



E(k, oS) dkdio 5i 



or just 

D i4 = D T S. 



T Oij 



(12) 



(13) 



with the eddy diffusivity Dt — w 2 t covi formed here only with 
the vertical turbulence intensity. 

More structure results if the turbulence is subject to a basic 
rotation. The expressions are given here only in the so-called 
r-approximation which is very close to the mixing-length ap- 
proximation (Kitchatinov 1986). We find 

Dij = D T (* x 5 y + * 2 0?fi°), (14) 

with f2° = as the vector parallel to the rotation axis. The 
components of the tensor are 

D rr = D T (*i +f 2 cos 2 60, (15) 
Dee = D T (*i + * 2 sin 2 6), (16) 
D r0 = D T *2sin6»cos6» (17) 

(see Hathaway 1984). The ^ ensures a latitudinal transport 
of the chemical composition. Without rotation $1 = 1 and 
$2 = 0. The rotational quenching functions result as 



(8) *i 



3 (_ , „ 2 wX arctanr2* 



1 + (fT 2 - 1)- 



(18) 
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and 

*2 



3 + (ft* 2 + 3) 



arctan ft* 



(19) 



4si* 2 v " ' v " ' ft* 

with the Coriolis number ft* = 2r corr O. 

2.2. Horizontal turbulence 

As usual we consider the temperature profile below the con- 
vection zone as stable so that any vertical fluctuations there 
are strongly suppressed. Any turbulence which possibly exists 
in the tachocline must be anisotropic. As a model of such an 
anisotropic turbulence, in this section a strictly horizontal tur- 
bulence is considered, i.e., there are, for zero rotation, no verti- 
cal fluctuations. 

For horizontal turbulence without rotation the diffusion ten- 
sor becomes 



Dk 2 



j 2 + D 2 k 4 



E 2 (k,uj) dkdui (Sij - g t gj) (20) 



without any vertical components. But with rotation and within 
the r-approximation we find 



D; 



= sD T j<% l*3 + 3(g ■n°y^ 2 ) :-!*_>(>;. 07 



(21) 



which yields in spherical coordinates the tensor components 

D rr = sD T (*3 + *4 + cos 2 61(6*2 + 2*5)) , 



D ee = sD T (* 3 +3* 2 ), 
D rB = D Br = -sD T sin cos 0(3*2 + *s)- 



(22) 



The amplitude s is introduced as the intensity of the horizontal 
of the intensity of the isotropic turbulence 



turbulence in units 

in the bulk of the convection zone 



*3 = 
*4 = 
*5 = 



ne duik 01 tne convection zone. 
The rotational quenching functions result as 

arctan ft* 



- (ft* : 

3ft* 2 



4ft* 2 
3 

Ail* 2 
3 

4ft* 2 \n* 2 + 1 



•7)- 



ft* 2 - 
ft* 2 + 3 



ft* 

(ft* 2 



+ (n* 2 - 3) 



. arctan fT 
arctan Q* 



(23) 



Without rotation * 3 = -* 4 = 1 and * 5 = 0, only the 
latitudinal diffusion coefficient Dee exists. With rotation, the 
appearance of the vertical diffusion D rr and the off-diagonal 
components D r e — De r can be observed (Fig. ||), which are 
playing - as we shall see - key roles in the mean-field equation 
for the large-scale concentration distributions in rotating turbu- 
lence fields. One can demonstrate that this equation for all f2* 
remains of the elliptical type. 

Note that the influence of rotation is so strong that for fT ~ 
1 the vertical diffusion - which is only rotationally created - 
equals the latitudinal diffusion. Note also that the new diffusion 
coefficients are strongly ^-dependent, resulting in ^-dependent 
profiles of the mean concentration (see below). 




4 6 
CORIOLIS NUMBER 



10 



Fig. 2. The diffusion tensor components (|22|) for midlatitudes 
(0 = 45) in their dependence on the Coriolis number fT. Note 
that for fast rotation (ft* > 1) all components are rotationally 
quenched 



3. The model design 

The mean-field diffusion equation without large-scale circula- 
tion reads 



dC_ 
dt 



)_d_ 

r 2 dr 



D, 



dC_ 

dr 



1 







r 2 sin 9 89 
The boundary conditions are 

C = 



D r 


9 <9C 


)} 


r 


89 




dC 

ae + 


rD r6 


dC 
dr 



(24) 



(25) 



at the lower boundary (xq = 0.6), where the lithium may be 
destroyed, and 



Drr dC + D ll dC =0 
dr r 89 



(26) 



at the solar surface (x = 1) where no radial flux may be al- 
lowed. 

For the correlation time in the definition of f2* a radial pro- 
file is used very similar to the radial profile of the turnover time 
in the mixing-length theory of the solar convection zone, i.e. 



n* = n 



7/6 / 1 



l-Si 



3/2 



(27) 



where ft* is the ft*-value at Xi = 0.7. Within the surface layer, 
x > 0.95, the Coriolis number ft* was put to zero. 

Eq. ( p4| ) is used in dimensionless form, distances measured 
in units of radius, R, and times measured in units of the dif- 
fusion time, Tdifi = R 2 /Dt- We start with a model close to 
that of Vincent etal. (1996). The diffusivity Dt may be con- 
stant in the whole integration region, xq < x < 1. Between 
the burning zone at xq and the bottom of the convection zone, 
x;, there is a turbulence field with motions only in the horizon- 
tal directions (Fig. ||). Its intensity is given by the parameter s. 
Following Spiegel & Zahn (1992) we shall refer to the region 
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Fig. 3. The overall turbulence model. Below the convection 
zone with isotropic turbulence there is a stable 'tachocline' 
layer with a horizontal turbulence field. At its bottom lithium 
is burned 



Fig. 5. The decay of chemicals at the surface of the convection 
zone for isotropic turbulence in the whole computational re- 
gion, xq < x < 1. No rotation (solid), with rotation (fi* = 6, 
dashed), extra horizontal turbulence in the tachocline beneath 
the convection zone (il* — 6, s = 100, dotted) 





Fig. 4. Isolines of the chemicals without (LEFT) and with 
(RIGHT) rotation feedback on the turbulence, 57* = 6. In both 
cases s = 100. Note that the rotation produces latitudinal pro- 
files in the composition 



xq < x < Xi as the solar 'tachocline'. The probing of its tur- 
bulence with respect to the lithium problem is the main scope 
of the present study. 

The general influence of rotation on the spatial distribution 
of chemicals inside the star is given in Fig. [i]for both ft* = 
and f2* = 6. The rotational influence produces a distinct de- 
pendence of the concentration on latitude. Figure || shows the 
decay of chemicals at the surface. The solid line corresponds to 
non-rotating turbulence for any value of s. The same result is 
obtained even for s = 0, i.e. without horizontal turbulence. The 
result of Michaud & Charbonneau (1991) and Charbonneau 
(1992) is confirmed in that 'horizontal turbulence cannot re- 
duce the effect of a given vertical turbulent diffusion coeffi- 
cient'. 

The dashed curve describes the depletion of the chemicals 
under the action of isotropic but rotating turbulence. Note the 
rotational slow-down of the chemical depletion. Faster deple- 
tion, however, is provided by the extra influence of rotating 
anisotropic turbulence beneath the convection zone. By this re- 
sult the influence of the rotation on the turbulence field is re- 
flected. The basic phenomenon is the formation of a latitudinal 
profile of the concentration as shown in Fig. f|. 



We have to take from Fig. |]that the additional anisotropy 
due to the basic rotation acts in an unexpected way: rotation 
suppresses the mixing of isotropic turbulence but, on the other 
hand, it enables a horizontal turbulence (which without rotation 
would not be active) to accelerate the vertical mixing. The rota- 
tion in the mixing process which is demonstrated in Fig. || has 
a very complex character: it both suppresses and deforms the 
turbulence. The interplay between both procedures yields the 
resulting effect. We shall see that for slow rotation enhanced 
diffusion dominates while for fast rotation the diffusion reduces 
the suppression of turbulence. As a result for a given rotation 
rate a maximal diffusion rate exists between slow and fast ro- 
tation. 

4. Lithium depletion and the dynamo problem 

Now the turbulence in the convection zone is considered as 
isotropic turbulence under the influence of rotation described 
in Section 2.1. Below the convection zone the turbulence may 
be so anisotropic that vertical motions do not exist. The ba- 
sic rotation, however, produces off-diagonal components in the 
diffusivity tensor as described in Section 2.2. Thus we have 2 
free parameters for the tachocline turbulence: intensity s of the 

(2) 

horizontal motion and correlation time r CO r r of the eddies. For 

(2) 

s — > the turbulence completely disappears and for T CO rr — > 
the rotational influence disappears. In both limits the decay 
time of the chemicals at the solar surface must become infinite. 

More information is given in Fig. |6j where the radial depen- 
dencies of the concentration profiles inside the star are given for 
three different latitudes and two different horizontal intensities. 
The diffusion is now much slower than with isotropic turbu- 
lence in the whole computational, space as in Fig. |5|. The differ- 
ential concentrations are produced in the tachocline and are di- 
rectly imprinted at the surface and are therefore not screened by 
the convective pattern within the convection zone. This is the 
same situation as for differential rotation but different from tur- 
bulent heat transport (cf. Stix 1981). The time-dependence of 
the surface differences of the concentration is so weak that our 
snapshots are always characteristic the situation. In Fig. ^ the 
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Fig. 6. Radial profiles of the chemical concentration along the 
equator (solid), mid-latitudes (dashed) and the poles (dotted) 
for s = 1 (LEFT) and s = 0.01 (RIGHT). Only for small 
horizontal intensities does the pole-equator difference become 
rather small. The time-dependence is now very weak, fi* = 6 
at the base of the convection zone and also within the tachocline 
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Fig. 7. Depletion time in units of the diffusion time. For weak 
horizontal turbulence intensities it approximately varies with 
1/s.Cl* =6 at the base of the convection zone and also within 
the tachocline 

decay time of an initially uniform concentration is plotted for 
various horizontal turbulence intensities in the tachocline. For 
s = 0, of course, there is no depletion, but the depletion time is 
strikingly short for nearly homogeneous turbulence fields with 
(say) s = 1. The correlation time profile and the basic rotation 
are taken from a solar model by Stix & Skaley with Q* = 6 at 
the base of the convection zone. We find that rotation is highly 
effective in transporting the chemicals to the burning zone at 
Xq. A scaling such as 



Tdc 



Tdiff/s 



(28) 



is found compatible with the idea that the diffusion of hori- 
zontal turbulence under the influenc of rotation with Q* ~ 1 
can be considered approximately as diffusion of isotropic tur- 
bulence. In reality, both the characteristic times differ by 7 or- 
ders of magnitude. Therefore, the horizontal rms velocity of the 
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Fig. 8. Depletion time Td oc (in units of the diffusion time) for 
horizontal turbulence intensities with s = 1 strongly depends 
on the correlation time of the tachocline turbulence, taken here 
for a rotation period of 25 days 

tachocline turbulence must not exceed 10 - ( 3 — 4 ) of the con- 
vection zone turbulence. That there is still lithium at the so- 
lar surface is only compatible here with the given concept of 
rotating turbulence with very small s of the tachocline turbu- 
lence. Consequently, the horizontal turbulence intensities must 
be very low, i.e. of the order of 1 cm/s. There is certainly no 
chance of maintaining magnetic fields of the order of kGauss 
by such a slow flow. 

But there is also the possibility that the tachocline turbu- 
lence has a rather short correlation time, e.g. similar to that of 
granulation at the top of the convection zone. Then the Coriolis 
number, f2*, is smaller than unity and the rotational influence 
can only be very small. For a given intensity (s = 1) the re- 
sulting depletion times are shown in Fig. [| We find a rela- 
tion Tdec/Vdig ~ 10 9 /r^ orr with T corr in seconds. Similarly, for 
s = 0.01 we find 



10 1 



T de 



"Tdiff 



(29) 



from our simulation. In the latter case with r corr ~ 10 3 ' s the 
lithium decay time will approach the order of 10 9 yrs. The cor- 
responding Coriolis numberis always smaller than 0.05. Again, 
with such a small value there is no hope of forming an effective 
a-effect for an appropriate dynamo process. In this approach, 
we find a close relation between the lithium problem and the 
theory of the solar dynamo. Within this framework, the obser- 
vations do not favour the existence of an overshoot dynamo 
below the bottom of the convection zone. 

5. Rotation rate variations 

For a given stellar model the basic rotation also forms a free 
parameter which can be varied. One can vary the basic rota- 
tion rate such as it varies in a young stellar cluster for all stars 
with the same age. The question is how the lithium varies. The 
same turbulence model can be used within the sample. We take 
s = 0.01 and TcoL =12 days. Figure plshows the results. Due to 
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Fig. 9. The general situation for young stellar clusters: The 
same turbulence model with s = 0.01 but the angular veloc- 
ity varied. The surface concentration at the equator is shown 
after 100 diffusion times. Due to the rotational quenching of 
the eddy diffusivity, the fast rotators possess the highest sur- 
face concentration of chemicals 

the rotational quenching of the eddy diffusivity (see Boubnov 
& Golitsin 1995), the faster the rotation the more lithium re- 
mains at the surface. The rotational quenching of the eddy dif- 
fusivity dominates the effect of the large latitudinal gradients of 
the concentration. The slower rotators are thus more effective 
in mixing the chemicals downwards. The plot, however, does 
not display that for still slower rotation, in the limit for fl* < 1, 
the mixing again becomes slower and slower and therefore the 
decay of the primordial lithium, too. The observation of the 
lithium-rotation correlation in young stellar clusters seems to 
confirm our prediction. 

6. Conclusions 

Rotation easily produces off-diagonal elements and other new 
terms in the diffusion tensor. We have shown with our com- 
putations that under the rotational influence, even a horizontal 
turbulence field will be very effective in transporting chemical 
components in the vertical direction. Two main consequences 
of this finding have been formulated, i.e. 

- even a strictly horizontal turbulence in the solar tachocline 
must be very weak (or very short-lived) to ensure the ob- 
served slow diffusion of lithium downward to the 2.6 Mio 
K layer where it burns, 

- for G-stars in young stellar clusters a correlation should ex- 
ist between the lithium abundance and the actual rotation 
rate. The faster the rotation, the stronger the eddy diffusiv- 
ity is quenched and the higher the lithium concentration in 
the convection zone must be. 

Our concept also provides statements for the long-term evo- 
lution of the lithium abundance under the influence of the gen- 
eral spin-down of solar-type stars. To this end a time-dependent 
code (which is in preparation) must simultaneously include 
both the time and depth-dependent effects of lithium burning 
and the decay law of the rotation rate with time. 
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